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C^ ' Abstract 

We propose a framework for calculating two-loop Feynman diagrams which 
appear within a renormalizable theory in the general mass case and at finite 
external momenta. Our approach is a combination of analytical results and of 
high accuracy numerical integration, similar to a method proposed previously 
[10 1 for treating diagrams without numerators. We reduce all possible tensor 



structures to a small set of scalar integrals, for which we provide integral rep- 
resentations in terms of four basic functions. The algebraic part is suitable for 
implementing in a computer program for the automatic generation and evalu- 
ation of Feynman graphs. The numerical part is essentially the same as in the 
case of Feynman diagrams without numerators. 
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Abstract 

We propose a framework for calculating two-loop Feynman diagrams which 
appear within a renormalizable theory in the general mass case and at finite 
external momenta. Our approach is a combination of analytical results and of 
high accuracy numerical integration, similar to a method proposed previously 



[10 1 for treating diagrams without numerators. We reduce all possible tensor 
structures to a small set of scalar integrals, for which we provide integral rep- 
resentations in terms of four basic functions. The algebraic part is suitable for 
implementing in a computer program for the automatic generation and evalu- 
ation of Feynman graphs. The numerical part is essentially the same as in the 
case of Feynman diagrams without numerators. 



1 Introduction 

The importance of developing techniques to calculate radiative corrections can hardly 
be overemphasized. The measurements at existing high energy colliders, like LEP 
and the SLC, are clearly sensitive to one-loop electroweak effects. As their precision 
increases, analyses at two-loop level may become necessary [Q. 

The problem of calculating one-loop Feynman diagrams was in principle solved 
completely long time ago by 't Hooft and Veltman [0], and by Passarino and Veltman 
1^. For a review, see for instance ref. |Q, and refs. |^, Q for some applications. As it 
is well-known, all one-loop diagrams are expressible in terms of Spence functions. In 
the general mass case, however, the numerical evaluation of the resulting expressions 
may be tricky |0]. This is because one needs to take care that the functions involved 
remain on the correct Riemann sheet, and because one has to control potentially large 
numerical cancellations. 

1 



No similar general solution was available so far at two-loop level in the general 
massive case, although a lot of effort was devoted to solving two-loop massive integrals 
at finite external momenta. There are two issues involved, given a general diagram 
with particles with spins: First, the numerator of the amplitude usually has non- 
trivial tensorial structures. It is convenient for us to decompose the amplitude into 
various tensorial constructs made of external momenta and metric tensors, multiplied 
by scalar functions, which in turn must be related back to the original Feynman 
integral by certain projections. Second, one should advance an algorithm by which 
one can evaluate these generalized scalar integrals, which have non-trivial numerators, 
in a systematic and, better yet, universal way. 

Regarding the evaluation of two-loop scalar integrals, a lot of work was done in 
this direction lately, and it has become clear that in general the integrals involved 
at two-loop level are not expressible analytically in terms of well-known and easy to 
evaluate functions like polylogarithms, as is the case with one-loop diagrams. Partly 
numerical approaches seem unavoidable in the general mass case. Only in some rare 
situations is an analytical solution in terms of special functions known. Such an 
example is ref. ||^, where a certain two-loop self-energy diagram was related to 
the Lauricella functions. A large number of approaches were proposed which work 
for specific topologies or mass combinations - see, for instance, refs. — [^ and 



references therein. In most cases, these works deal with Feynman integrals with 
simplified couplings and numerators. Some of these approaches can yield precise 
numerical results for the diagrams to which they apply; by combining several such 
methods, one was able to apply them to actual physical processes. While they do hold 
promise to be extendable to cover more complicated scalar integrals which may appear 
in processes involving particles with spin, derivative couplings, and more external 
momenta, it is fair to say that this program has not been completed so far. Moreover, 
for calculating a physical process, such approaches would involve separate methods 
to treat individual graphs. This implies a considerable amount of work which grows 
very fast with the complexity of the process considered. 

Fewer results exist concerning the reduction of two-loop Feynman graphs to scalar 



integrals. In ref. |2^ it was shown that the two-loop self-energy diagrams can be 
reduced to scalar integrals without numerators. This result applies only to two-point 
functions. It was used in conjunction with numerical integration in ref. |p4[ for 



calculating certain two-loop contributions to the p parameter. The approach of ref. 



[23 1 only refers to propagator diagrams. Ref. EH], which also deals with propagator 



diagrams, might be extendable to some more complicated diagrams, like certain 3- 



point functions in simplified kinematic cases, and the recurrence relations of ref. [^ 
perhaps can be solved numerically in more complicated cases, as suggested in the 
Conclusions of ref. ||2^ . 



Compared to the approaches briefly discussed above, in this paper we present a 
general universal framework which applies to any n-particle 2-loop diagram with arbi- 
trary renormalizable dynamics and any combination of masses and external momenta. 



Strictly speaking, our formulae can be used directly for calculating one given diagram 
only if that diagram is free of mass singularities. However, in actual calculations this 
can be often circumvented by isolating the singularity analytically at the level of the 
integrand, by considering infrared finite sums of diagrams, or by introducing a mass 



regulator, as was possible for instance in ref. [|I^. In our approach, all diagrams are 
treated using the same algorithm, and for this reason the whole calculation, including 
the generation of the relevant diagrams, can be authomatized. This was done in ref. 



10 1 in a simplified case of if — > WW decay, which involves no tensor structures in 



the numerator, and is a subset of what we shall discuss in this paper. This can be 



compared with ref. [T^, whose authors used separate methods for individual dia- 
grams to perform the same calculation, and which is so far the state-of-the-art of 
physical calculations by using these methods. 

Within the framework which we propose, we shall show that in any renormalizable 
theory all two-loop diagrams can be expressed in terms of a few scalar functions. To 
our best knowledge, these functions cannot be expressed analytically in terms of 
known functions, in a way which would facilitate their numerical evaluation. There 
should be in principle a connection with generalized hypergeometric functions, but it's 
not clear that this would lead to an efficient evaluation of these functions. Instead of 
looking for an analytical result, we derive one-dimensional integral representations for 
these scalar functions. They can be expressed in terms of four simple functions. The 
structure of these integral representations is a generalization of the results derived 



in ref. [10 for the case of two-loop diagrams without numerators. Such integral 



representations can be used for calculating numerically the necessary functions very 



efficiently. This was shown in refs. |jTO|, |T^ for two-point functions and in ref. [Q 
for three-point functions. 

Compared to the case of diagrams without numerators, two new but similar struc- 
tures may appear in a renormalizable theory, like the standard model. We note that 
in nonrenormalizable theories more structures are allowed. As it will become clear 
from our discussion, these additional structures can be treated along the same lines, 
and will result in similar functions. 

Added Note After the submission of this article, a related work by O.V. Tarasov 
appeared [^. This work makes use of some special properties (dimensional shifts), 
due to Schwinger representation of propagators and was able to express any two-loop 
propagator diagrams by a minimal set of four master integrals. It is further asserteded 
that this result does not depend on the renormalizability of the theory. On the other 
hand, our assertion is that any two-loop n-point diagram can be expressed by four 
basic functions, as long as the theory is renormalizable. We have not specifically 
looked into propagator diagrams to see if there are extra symmetries or recurrence 
relations such that this result still holds for nonrenormalizable theories. Furthermore, 
we were not primarily interested in identifying a minimal set of basic functions, but 
rather in obtaining a framework which leads to expressions convenient to evaluate 
numerically. Thus, our remark that there is no finite set of basis integrals for non- 



renormalizable theories should be taken only as a most direct inference from our work 
so far. 



2 Reduction to scalar invariants 

In this section we show that the calculation of two-loop Feynman diagrams with any 
tensor structure at the numerator can be reduced to the evaluation of a class of scalar 
invariants. 

Any two-loop diagram without numerators can be written as an integral over 
scalar integrals of the type ||T0| : 



(p2 + mf )"i (g2 + r7ii)°2 [(r + kf + mjj^'^ ' ^ ^ 

where all momenta are Euclidian. 

This form is obtained after combining all propagators which contain the same loop 
momentum p, g or r = p + g by introducing Feynman parameters. At their turn, all 
G functions can be derived from two basic functions JF and Q. For a discussion of the 



properties of these functions see refs. |]10|, |Ij]. To calculate a Feynman diagram, one 
typically has to perform a numerical integration of. eq. 1. The maximum dimension 
of this integration is the number of propagators of the Feynman diagram minus three. 
A Feynman diagram may contain non-trivial tensor structures at the numerator. 
For all these cases, tensor structures in the original Feynman diagram will result in 
tensor integrals of the type: 

^ "^ (p2 + mf)"i (g2 + rni)"2 [(r + k^ + mi]"3 ' ^ '' 

The calculation of such tensor integrals can be reduced to the calculation of scalar 
integrals of the following form: 

rf> dn, ip-kn<i-kf 

(p2 + mf )°i (g2 + mi)°2 [(r + k)^ + m|]°3 ' ^ ^ 

To show this, let us introduce the transverse components of the loop momenta p 
and q: 



P± - P ^2 '^ 

ql = q'- ^fc" • (4) 



With these notations, eq. 2 resuhs in a sum of terms with the following structure: 

J {p^ + mlY^ (g2 + mlY^ [{r + ky + mlY^ ^ ' 

where P is some poynomial function oi p ■ k and q ■ k. One can convince oneself 
that the l/k"^ singularities introduced in these expressions by the decomposition of 
eqns. 4 are superfluous. They disappear completely when one adds the transverse 
and longitudinal parts together. We shall show this explicitly in future publications 
devoted to calculations of specific processes. 

Integrals of this type vanish for odd j'. This is because the tensor /^^ '^j' is 
transverse with respect to all its indices. At the same time, the only vector available 
for constructing fi'^'j' after the p and q integrations in eq. 4 are carried out is the 
external momentum k'^. If j' = 2n + 1, after using k and the metric tensor g^'^ to 
support 2n transverse indices, there is one free index left, which must be carried by k. 
However, this last index cannot be made transverse, and thus our assertion follows. 

As an example, one has the following relation: 

p^P{p-k,q-k) 
dp a q ■ 



{p^ + mf )"i (g2 + m|)"2 [(r + ky + m|]"3 



tl I d-p d-q p-kP{p-k,q-k) 

P y (p^ + mf )"i (g2 + m|)°2 [(r + ky + m§]"3 ^ ^ 

This is because the following relation holds: 

^ ^ (j92 + mf )"i (g2 + mlY^ [(r + kf + m^J^s ~ ' ^ ' 

If j' in eq. 5 is even, one can always decompose fi'^'j' into scalar invariants of 
the type in eq. 3 by using the transversality of I>^^"'^^J' . 

For instance, apart from a redefinition of the loop momenta p and g, two combi- 
nations are possible for j' = 2: 

{p^ + mf )°i (g2 + mlY^ [(r + kf + m§]"3 ^^ p '' 

dyd-q f^'iTnP-K'i-k) ^ (,...._ ^^!:^)5 (8) 



where 



^ = ^U-pd^q plPip-k,q-k) 

n-lJ {p^ + mJY^ (g2 + m|)°2 [(r + A;)^ + m|]"3 

1 I d^pdTq (P± • <1A-)P{P -k.q-k) 

n-1 J {p^ + mJY^ (g2 + m|)°2 [(r + ky + m|]"3 



5 



The scalar invariants A and B can be further reduced to functions of the form of 
eq. 3 by using the relations: 

p1 = P^ - -j^{p- kf 
P±-q± = p-q--^{p-k){q-k) 
p.q = -(r -p -q) 

and by partial fractioning the resulting expressions. 

The corresponding relations for j' = 4 are given in Appendix A. 

Obviously, this procedure can be extended to higher order tensor structures which 
may occur in more complicated calculations. At this point, any two-loop Feynman 
diagram can be expressed in terms of scalar integrals of the type of eq. 3. 



3 Recursion relations 

We have shown in the previous section that the problem of calculating any two-loop 
Feynman diagram reduces to treating scalar integrals of the type: 



rprn ip-kn<i-kr .,,. 

(p2 + mf )"i (g2 + mi)"2 [(r + k^ + mi]"3 ' ^ ^ 

where ai, 0:2, 0:3 > 1 and a,b > 0. 

Not all ^0^0203 ^^^ independent. There are recursion relations between functions 
with different indices. In this section we show that all V^^^^^,^ functions which can 
appear in a two-loop calculation within a renormalizable theory can be derived by 
differentiation from a limited number of basic functions. 

As a first step, one notices that one can increase the indices ai, a^ and aj, of 
^01^02 Q3 by one unit by differentiating with respect to the mass arguments: 

1 d 
Kr+ia2a3i'rni,m2,m3;k'^) = - — -^V^^^^^.^imi,m2,m3;P) , (11) 

and similarly for 0:2 and 0:3. 

Next, by differentiation with respect to the external momentum k^, one obtains 
the following relations: 



-pa+l b __ 

' ai+l«2a3 2 ry 



2*^A -(« + .) 



jyab _|_ "^ pg-lfe 

ai a2 "3 2rv "^ "^ °3 



v: 



ah+l 

ai a2+lcK3 



2^2 



2k' 



_d_ 
9F 



(a + 6) 






V. 



ab-l 



(12) 



and 



2t 



_d_ 



- (a + 6) 



v: 



ab 

ai 02 "3 



-2a, 



7 2pa6 

"" ' ai a2 13- 



>a+16 



>afe+l 



-1 ' ' ai 02 c«3+l aia2<^3+l 



(13) 

In eqns. 12, the functions P^i 02^03 °^ ^afa2a3' which may appear if either a or 6 
vanish, are defined to be zero. These relations allow one to increase the upper and 
the lower indices by one unit simultaneously. 

In the following we will call 01 + 02 + 0:3 — — 6 the "degree" of the function 
'Paia2a3- Equs. 11 cau be used for increasing the degree of V^^^^^^^. Eqns. 12 
essentially relate functions of the same degree, while increasing simultaneously the 
upper and the lower indices. 

Hence, it is important to look at the functions of lowest possible degree in order 
to identify a limited set of functions from which one can derive all the other. 

The allowed range of the degree of the functions P^f 0,2^3 which can appear in 
two-loop amplitudes depends on the theory under consideration. For instance, in the 
case of the linear sigma model the minimum degree is three M, O]. There is no 



upper limit, since the degree can be increased indefinitely by adding external legs to 
the diagram. 

In renormalizable theories, the degree of two-loop functions has always a lower 
bound. To see this, one notes that powers of the loop momenta appear in the numer- 
ator of a Feynman diagram only from the fermion propagators and from derivative 
couplings. If one demands the theory to be renormalizable, one can only have fermion- 
boson couplings with no derivative, and trilinear boson couplings which have at most 
one derivative. One can convince oneself that the minimum degree of a two-loop 
function is one, and is attained by the vacuum diagram with two fermion propagators 
and one boson propagator shown in fig. 1. By adding external bosonic lines or pairs 
of fermionic lines to a diagram one can only increase the degree of the diagram. 

Therefore, by using eqns. 11 and 12 one can derive all functions V^^^^^^ which 
may appear in renormalizable theories from those with oi = 02 = 03 = 1 and 
a + 6 = 0, 1, 2. Nevertheless, we choose to use the following set instead: 



degree = 4 
degree = 3 
degree = 2 
degree = 1 



ppo 



211 



plO 



211 



p2 



211 



p3 



211 



pOl 



211 



pll 



211 



p21 



211 



pp2 



211 



pi 2 



211 



pps 



211 



(14) 



The reason is that for these functions one can obtain simpler integral representations 
than for the functions with Oi = 02 = 03 = 1. The functions with Oi = 02 = 03 = 1 




Figure 1: Two-loop vacuum diagram with one boson and two fermion propagators, 
which leads to a P function of minimum degree (one). Vacuum diagrams are expressed 
analytically in terms of Spence functions, but one can obtain nontrivial V functions 
of degree one by adding external bosonic lines to this diagram. 



are then derived from the set of eq. 14 by using the "partial p" operation of 't Hooft 
and Veltman 01. 



By applying partial p to V^^^^ ^^ 5 o^^ finds the following relation: 



v: 



ab 

ai 02 03 



n- {ai + a2 + as) + (a + b)/2 



{amlV:'^+ia,a, + «2m2p»f,^+i 



03 



+«3m2p^f^^^3^,-- 



2k' 



_d_ 

'die 



[a + h) 



' ai a2 "3 I ' y / 



and this can be used together with eqns. 12 or 13 for expressing the functions with 
tti = a2 = ^3 = 1 in terms of the set in eq. 14 when they are needed. 

Note that not all functions in eq. 14 are independent. Some of them are related 
by loop momentum redefinitions. For instance, the following relation holds: 



01 



V2ii{mi,m2,m^]k ) +V2ii{mi,m2,m^]k ) 



V2ii{mi,ms,m2;k ) + k V2ii{mi,m3,m2;k )^ . (16) 

At the same time, some of these functions may actually never appear directly from 
a Feynman diagram, as is the case with P^ii- However, they may appear indirectly 
when using eqns. 12, 13 or 15. For this reason we prefer not to discard them at this 
point. 

The set of functions in eq. 14 is an extension of the functions JF and Q of ref. ||1CI|| . 
!F and Q are enough for treating all diagrams without numerators. This is the case for 
instance with the linear sigma model. As already mentioned, the lowest degree which 
may appear in the linear sigma model is three. Therefore, only the first two lines of 
eq. 14 are involved. In fact, P211 — S, and P211 — ~^ ■ The remaining degree three 
function, P^ii, never appears directly in this case from a Feynman diagram. 

One also notices that eq. 15 is a generalization of eq. 11 of ref. [O]. 



Finally, a remark on nonrenormalizable theories is in order. One can construct 
two-loop Feynman diagrams of degree lower than one if one allows for nonrenormal- 
izable interactions. As an example, the two-loop sunset self-energy in a four-fermion 
interaction theory is of degree zero. In nonrenormalizable theories it is not always 
possible to restrict the degree of the V^^^^ ^^ functions and thus to identify a finite 
set of basic functions which generate all the other. In the following section we will 
derive one-dimensional integral representations for the functions of eq. 14. As it 
will become clear, it is possible to derive similar integral representations for functions 
with a + 6 > 3 as well. Thus, specific calculations in nonrenormalizable theories can 
be performed, too. One only needs to introduce additional functions along with those 
of eq. 14. 

4 Integral representations 

We have seen in the previous sections that any two-loop Feynman diagram in a 
renormalizable theory can be expressed in terms of the ten functions of eq. 14, not 
all of them independent. 

As said, our aim here is to express the ultraviolet finite parts of a calculation 
in a general kinematical situation in a form which would facilitate their numerical 
evaluation. The reader may know that, already at the level of two-point functions, 
the self-energy sunset diagram is known to be related to the Lauricella functions 0. 
These functions are not straightforward to evaluate numerically. The best way to 
do that appears to be by means of integral representations [|19|. For more compli- 



cated two-loop diagrams no analytical results are known so far. Besides, we want to 
deal with more realistic situations, where non-trivial numerators are present in the 
integrands. 

Hence, in this section we derive one-dimensional integral representations of the 
ten functions in eq. 14. Such integral representations can be used for a fast and 
accurate numerical evaluation of these functions. 

It is convenient to replace the set of ten functions 7^211(1711, 1712,1713; k"^) of eq. 14 
by the following equivalent set 7ii{mi, 7712,1713; k"^): 

1 



-pOO _ 

' 211 


^7i, = 


<rp<rq 


-pio 

' 211 


-^7i2 = 


j (TpcTq 


' 211 


-^3 = 


j <rp<rq 


p20 _ 
' 211 


-^7ii = 


[ <rp<rq 



[{p 


+ ky 


p ■ k 


mi) (r2 


+ 


^i) 


[{v 


+ ky 


q ■ k 


mi) (r2 


+ 


ml) 


\{p 


+ ky 


{vkf-h 


mi) (r^ 
^p^ 


+ 


ml) 



[{p + k)"^ + mf]2 (g2 _|_ TT^i) (r^ + m|) 



' 211 


-n, = 


/ t/>rf"g 


-p02 
' 211 


-^6 = 


j <rp<rq 


-p30 . 
' 211 


-^n, = 


j (TpcTq 


p21 , 

' 211 


-^8 = 


j <rp<rq 


pl2 
' 211 


-^9 = 


j <rp<rq 


-p0 3 
' 211 


-Hio = 


fdyrq 



(p-k){q-k) --k'^{q-p) 



[{p 


+ ky 


+ rr 


kf- 


+ mi)(r2 


+ mi) 


Up 


+ ky 
{p 


+ m 


3 

n+2 


+ mi)(r2 
k'^p^{p ■ k] 


+ mi) 


Up 


+ kY + m 
{p-k)\q 


'i? iq' 
k)- 


42^V(g 


+ mi) 
k) 



[{p + A;)2 + mi]2 (g2 + 77^2) (r^ + ml) 

2 ^P[2(p-g)(g-A;) + g2(p 



,^Jp-A;)(g-A;) ^_^^ 



[{p + A;)2 + mf]2 (g2 _)_ ^2^ (^7-2 _|_ ^2^ 

(g-^)^~42^^g^(g-^) 
[(p + ky + mf]2 (g2 + ■m'2) (r^ + mg) 



(17) 



For compactness, we omitted in the above formulae the mass and momentum argu- 
ments of the functions V2ii{mi,m2,m3;k'^) and 7ij(mi,m2,m3; /c^). 

We choose to consider this set of scalar integrals because their ultraviolet be- 
haviour is logarithmic. This is necessary for deriving simple integral representations 
of their ultraviolet finite parts. It is straightforward to express the P2 1 1 functions in 
terms of the Tii functions. The necessary conversion formulae are listed in Appendix 
B. 

Simple, one-dimensional representations for the functions ?^i and 7^2 were found 
in ref. [0. As it turns out, it is possible to derive similar integral representations for 
the other functions, H3 — Tiio, as well. We discuss some technical details of deriving 
such integral representations in Appendix C 

By using the methods discussed in Appendix C, one finds the following expressions 
for the functions Hi (e = n — 4): 



Hi 

^2 

Ha 



TC 






4;,2 



vr^fc 



2^1.1 . s^l3 7r2 7„^ 2 , 



n 



*e 



e'2 

-^ - -(t - 7mi) -777 + 777 ^ + ^ + ^3 



7r\ky 

7l\ky 



3 132V 
2e2 e 2 



8 12 
13 7r2 
16 "^24 4 
175 7r2 372 3, ■ 

1 ^ -Lull. + _/;, 

96 16 4 4 



137„, 175 



TT 



37 



2 

mi 



4e2 e 4 192 32 

1 1,1 7mi. 



3/, 



2e2 e^24 



19 , TT^ 7^, 7mi , 3 ' 
32^48 24 ^ 4 ^ 4 ' 



10 



n. 



n. 



Ho 



n 



10 






1 1 5 



Imi, 



287 7r2 57™, 7™i 1 



192 24 24 2 2' 

.2 cr^, ^,2 



/i7 



mi , 7™, ^ 1^^ 



1 1,5 7„, , 287 TT^ 57, 
2e2 eM8 2 ^ 384 48 48 4 ' 2 

^ _ 1 /^ Tm^x 95 TT^ 7^^ 7^^ 1 

3e2 6^4 3 ^ 192 72 24 6 2^ 

4? + 7^96 + ~^"768 + 96 + ^ + ~ + 2^^° 



'1^ 



Again, to simplify the notations, we omitted in the above formulae the mass and 
momentum arguments of the functions 7ij(mi, m2, ms; /c^) and hi^rrii, 1712, nis-jk'^). 
7mi = 7 + log {nm\/ n\) , where 7 = 0.577216 is the Euler constant and /xi is 
the 't Hooft mass. The ultraviolet finite parts hi{mi,m2,'m^;k'^) of the functions 
7ii{mi,m2,m^; k"^) have the following one-dimensional integral representations: 



/ii(mi,m2,m3; k" 



h2{mi,m2,ms;k' 
/i3(mi,m2,m3;/c^ 
/i4(mi,m2,m3; A;^ 
/i5(mi,m2,m3; A;^ 



hfi{mi,m2,mz]t 



/i7(mi,m2,m3; k' 

hs{mi,m2,ms;k'' 

/i9(mi,m2,m3;A;^ 

/iio(mi,m2,m3;A;^ 



dx g{x) 



dx [g{x 


) + /iW 


dx [g{x 


) + fi{x) 


dx [g{x 


) + fi{x) 


dx [g{x 


) + fi{x) 


dx [g{x 


) + fi{x) 


dx [g{x 


) + fi{x) 


dx [g{x 


) + fi{x) 


dx [g{x 


) + fi{x) 


dx [g{x 


) + fi{x) 



\{l-x) 

+ kx)] 

+ f,(x)]il-x) 
+ f,(x)]il-xf 

+ f2{x) + hix)] 

+ f2{x) + Mx)]{l-x) 

+ kx) + hix)]il-xf 
+ kx) + Mx)]{l-xf 



(19) 



All ten integral representations are built out of the following four basic functions: 



g{mi,m2,m3;k'^;x) = Sp{- ) + Sp{- ) + t/i log ^+ 1/2 log ^ 

1 - yi 1-^/2 yi-l 1/2-1 



/i(mi,m2,m3; A;^;x) 



l-/i2 



K^ 



+ yl log 



yi 
yi-i 



1/2 log 



1/2 
1/2-1 



11 



f2{mi,m2,m3;k ;x 



f3{mi, 1712, m3;k^;x) 



1 
4 






K^ 




+l/l log 7 + 2/2 log 7 

yi - 1 y^ - 1 



,(20) 



where we use the following notations: 



2/1,2 

A 



l + K^-fI^±VA 



2k^ 



'\^k' 



2\2 ,,(22 
/i ) + 4k /i 



AiK^r] 



(21) 



and 



/^ 



ax + 6(1 — s) 
x(l — a;) 

^2 ' — 2 5 



mt 



m 



2 



(22) 



In the above expressions, one special case must be treated separately, namely 
k^ = 0. One can convince oneself that in this case our approach reduces to the 
functions introduced by van der Bij and Veltman in ref. [Q. 

The functions g and /i were already introduced in ref. [0. They are sufficient 
for treating Feynman diagrams without numerators. The other two functions, /2 and 
/a, have a similar structure. They are needed for treating functions of degree two 
and one, respectively, which may appear for instance in a renormalizable theory with 
fermions. Clearly, it is possible to extend these formulae to functions of lower degree 
which may appear in nonrenormalizable theories. 

In order to calculate more complicated diagrams, the derivatives of the functions 
defined in eqns. 20 are needed. One may ask oneself whether this procedure does 
not introduce additional singularities in the integrands, for instance at x = 0, 1. 
One can convince oneself rather easily that this is not the case. This can be seen 
by noticing that near x = (the situation is the same at x = 1), the function e.g. 
g, behaves like logx^, which is integrable. It is easily seen that by differentiating 
with respect to one mass or external momentum the behaviour near x = will 
remain integrable. By differentiating the function an arbitrary number of times one 
cannot induce additional singularities in the integral representation. The derivatives 



12 



will remain integrable. Since there can be no fundamental difficulty related to the 
differentiation procedure, the problem of evaluating numerically the derivatives is 
merely a question of implementing the formulae in a computer program in a correct, 



numerically stable form. This was shown in previous works, e.g. refs. |llO| , p!2| , which 
involve a subset of the hi functions, and where the differentiation procedure was used. 
The relations 19 can be used for an efficient numerical evaluation of two-loop 
Feynman diagrams. We do not enter here into details regarding the numerical in- 
tegration. Let us just note that this was done in refs. [1^ — [12] for hi and /i2 in 



the case of two- and three-point functions, and the results agree with independent 
calculations |jl3[ — |jl5[. The other h functions have analytical structures similar to hi 
and /i2, so that the same numerical algorithms can be used to evaluate them. 

5 Conclusions 

We presented a method for treating all two-loop diagrams which may appear in a 
renormalizable theory in a systematic way. Formulae were given which allow one to 
express any two-loop diagram in terms of a small number of scalar functions. We de- 
rived one-dimensional representations of these scalar functions. They are constructed 
out of four basic functions which have quite simple expressions. These integral rep- 
resentations can be used further for an efficient evaluation of Feynman diagrams. 

A useful feature of our approach is that it is suitable for the automatization of 
two-loop computations. All algebraic manipulations needed to reduce a Feynman 
diagram to scalar invariants and to express it in terms of the four basic functions 
which we introduced are algorithmical. They can be encoded in a computer algebra 
program to treat automatically two-loop Feynman diagrams. 

The calculation of a Feynman diagram implies a numerical integration over hi 
functions or their derivatives. The maximum dimension of this final integration is the 
number of propagators minus three. Where derivatives of hi functions are needed, 
they are straightforward to obtain, for instance by using a computer algebra program. 

The numerical integrations, the details of which we do not discuss in this paper, are 
standard. Let us just note that the analytical structure of the functions involved in the 
integral representations is known [0. This allows one to identify their singularities, 
and to define a complex integration path by means of spline functions. Along such 
a path one can use an adaptative deterministic algorithm for a fast and accurate 
numerical integration |p!0| , |I2|. 

In this paper we provide only the formulae which are encountered in renormalizable 
theories. If calculations in nonrenormalizable theories are addressed, new structures 
may appear. However, they can be dealt with along the same lines, and lead to 
similar functions. The main difference is that in nonrenormalizable theories it may 
be impossible to isolate a finite number of two-loop functions from which all other 
can be derived by differentiation. In this case one would need to consider a specific 
process and to identify the necessary functions. 
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Appendix A 

We give here the formulae for decomposing a tensor integral I'^^"''^i' of the type in eq. 
5 for the case j' = 4 into scalar invariants. 

Apart from a redefinition of the loop momenta p and q, there are three possible 
cases if j' = 4: 

^n ^n p^Ip^^p^Ip^^P{p ■k,q-k) ^ 

^ ^ (p2 + mf )"i (g2 + m|)"2 [(r + ky + m|]°a 

(J J ( _y^2^i^2„M3M4 _|_ gMlM2y^M3y^M4 _|_ gM3M4y^Mly^M2 J^fJ-1 J^fJ-2 J^tJ-3 J^fJ-4, 

+ [(/Ul/i2/U3At4) ^ (/il/U3/i2/U4)] + [(/il/U2/i3/i4) ^ (1^11^41^21^3)] } (23) 



^n ^n g± P± P± P± -P(P ■k,q-k) ^ 

^ ^ (p2 + mf )°i (g2 + m|)°2 [(r + A;)^ + m|]"3 

1 



D I f -A;2^^1^2^^3/.4 ^ ^WM2p3p4 ^ ^M3M4^mp2 _ pi^M2^/.3^M4 

+ [(AtlAt2At3/U4) ^ {f^lfJ'3f^2fJ'4)] + [(/Ul/i2/i3/i4) ^ {fJ'lf^4f^2f^3)] } (24) 



^„ ^„ q^^q^^p^lp^^P{p-k,q-k) 



ijp' + mf )°i (g2 + m|)°2 [(r + A;)^ + m§]"» 
1^ (-Pg^'^^g^'-'^'^ + g^^'^^k^'k^"^ + g^'^^^k^'k"^ - ^k^'k^'k^^'k"^ 

+ -^ {[(AH/U2/U3At4) ^ (Atl/U3/i2/U4)] + [(/il/i2At3At4) ^ (/Ul/i4/U2/i3)] }} , (25) 



where 



C^ = - ., ^ .,J rprq- (Plfnvk,q-k) 



{v? - 1)A;2 7 (p2 + mf )"i (g2 + mi)°2 [(r + A;)^ + mi]"3 

I^ = ---i--/rf>rf"g (l>.-g.)l>iP(p-A;,g.A;) 



V? - 1)A;2 7 (p^ + mf )"i (g2 + mi)"2 [(r + kf + mi]"^ 
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E 



(n2-l)(n-2)A;2 



/ 



d>rf"g 



[nplql - 2(p_L ■ q±f] P{p -k^q-k) 
(p2 + mf )"i (g2 + mi)"2 [(r + A;)^ + mi]°3 



(n2-l)(n-2)P 



^„ ^„ [pW± -{n- l){p± ■ q±)^] P{p ■k,q-k) 
(p2 + mf )"i (g2 + m|)"2 [(r + A;)^ + m|]"3 



(26) 



The scalar invariants C, D, E and F can be further reduced to functions of the 
form of eq. 3 by using the relations: 



p1 



P±-Q± 
p-q 



P^ - -j^ip- ^f 



1 



P-<l-j^{p-k){q-k) 



(2 2 2\ 

[r — p — q ) 



and by partial fractioning the resulting expressions. 

Appendix B 

In section 4 we replaced the set of functions 7^2 ii ' a + & = 0, 1, 2, 3 by the functions 
Til — TYio- The functions Tii are free of quadratic divergencies, and for this reason 
they have simpler integral representations. The relations between 7^211 ^"^^ '^i ^^^ 
easily obtained by partial fractioning the eqns. 17: 



-pOO 
' 211 

-pio 

' 211 

-poi 

' 211 



' 211 



-pii 

' 211 



2 



-pu/ 
' 211 



3 



' 211 



Hi 

-n2 - k^Hi 

e 

n 



n^ + -{[{n- 1)P - ml]ni + 2(n - 1)7^2 + n°iM 

n^ + k^Hs +'^{{ml + ml-ml + k^)7ii + 27^2 - n°i°i 

+T2K)[ri(m2)-Ti(m2)]} 



Tie, 



n 



-mlHi + T2imj)Ti{ml] 



-n. 



3P 

n + 2 



n — 1 



r - mf k'Hi + [{n - l)r - mi]n2 + n7^4 - T^i 1 1 
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p21 

' 211 



-pi 2 

' 211 



-pUd 
' 211 



'^^2 r9 



n — 1 



01 



A;^ - mn ^3 - ^l^^ 



n — 1 

n(n + 2) 



ik 



2\2 



X 



{(m^ + m^ - m^ + e)ni + 27^2 - n°?i + T2iml)[Ti{ml) - Ti{ml)]} 



— TYq — k Tie 



n + 2 



,2k 



2^5 + ( m^)H2 + (mi + mi-ml + e)n^ + PJ'i^ 



n 



{k 



2\2 



n(n + 2) 



-Ti; 
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n + 2 



? - (n + l)m2 - m^ + A;^)?^! - ^0°, 

+T2K)[Ti(m2) -{n- l)Tr{ml)]} 
-mlHs + eT2{ml)Ti{ml] 



where Ti and T2 are the Euchdian one-loop tadpole integrals: 



(27) 



Ti(m' 



T2(m' 






p2 _|_ ^2 



— vr^ f vrm^j ^ F 



TT^ iTrm^ ) ^ r 



e\ 2m 



2/ 2 + e 

"2. 



(2^ 



(j9^ + m^)^ 

For simplifying the notation, we omitted in the above formulae the mass and mo- 
mentum arguments of the functions 7ij(mi, m2, m^; k"^) and 'P^^^^^.^{mi,m2, m^; k"^), 
and understand that these arguments appear in this order in all relations. 

The functions with ai = a2 = a^ = 1 which appear in the relations 27 can be 
calculated by partial p (eq. 15): 



V^i^{mi,m2,m3;k 



Pm(mi,m2,m3;A; 



n — 



^^vr,,-vll,-ml{n2 + en,] 



(mi,m2,m3; A;^ 



-7712^3(7722, mi, ms; k'^) + ml(H2 + H3) {1713,1712, mi; A;^)| 

f [(771^ + k'^)ni + 7^2] {rni,m2, m^; k"^) 

+m\'Hi{m2,mi,m3] k"^) + ml'Hi{m3,mi,m2] k'^)\ (29) 



Appendix C 

Here we make some comments on the derivation of the integral representations of 
the functions Tij. The derivation proceeds in the same way for all ten functions, and 
works for functions of lower degree as well. 

Let us consider for instance the function Tip.: 



H6{mi,m2,m3;k'^ 



cTpcTq 



{q-kf-lk\ 



ii.2^2 



[{p + kY + 7rif]2 (g2 _|_ 77^2^ (7-2 _|_ 77^2^ 



(30) 
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We introduce two Feynman parameters x and y, and integrate out the loop mo- 
menta p and q. Note that the quadratic divergence which is present in the P211 
integral cancels out in Tig. This is necessary for obtaining a well-defined integral 
representation of the ultraviolet finite part. Without this cancellation, part of the 
ultraviolet divergence of the diagram is transferred from the radial loop momenta in- 
tegration to the Feynman parameters integration; one is then left with a nonintegrable 
singularity in the x integration. 

After integrating out the q and p loop momenta, one finds (n = 4 + e): 



H6{mi,m2,m3;k^ 



A/u2\2-nl ^^Y(c^ " 



r 3 



TT^kY^ 



X / dx[x{l — x)]2{l — xY 
Jo 



Txm^ 



X / dy 
Jo 



i?2 



1-y 



-K-4 



(31) 



Here we introduced the following notation in addition to eqns. 21 and 22: 



L 
R 



2 2 
y K 



/i 



y{l-y)K'^ + y+{l-y)fj,^ = K^{yi ~ y){y - y2) 



(32) 



Eq. 31 displays already the structure of the function Tie- The first line of eq. 31 
contains the ultraviolet singularity, while the y and x integrals are finite. Hence, the 
X and y integrations need to be performed up to O(e^). Therefore, one expands the 
integrand in powers of e up to order 0{e^) and carries out the y integration. The 
calculation can be simplified considerably by using the relation: 

L y dR , , 

Finally, one obtains after the y integration the results given in eqns. 18 — 20. 
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